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Abstract 

Stochastic growth phenomena on curved interfaces are studied by means of stochastic partial 
differential equations. These are derived as counterparts of linear planar equations on a curved 
geometry after a reparametrization invariance principle has been applied. We examine differences 
and similarities with the classical planar equations. Some characteristic features are the loss of 
correlation through time and a particular behaviour of the average fluctuations. Dependence on 
the metric is also explored. The diffusive model that propagates correlations ballistically in the 
planar situation is particularly interesting, as this propagation becomes nonuniversal in the new 
regime. 

PACS numbers: 68.35.Ct, 02.40.Ky, 05.40.-a, 68.35.Fx 



1 



I. INTRODUCTION 



Growing interfaces appear everywhere in the natural world. Some of them are as impor- 
tant as the technologically motivated thin film deposition, medically relevant as bacterial 
colonies development, and physically interesting as fluid flow in porous media. Together with 
these, one finds a large variety of situations where the main driving mechanism is a growing 
fluctuating interface. Despite their diverse origin, all these phenomena have been studied 
within the common formalism of scaling analysis This theory classifies different growth 
phenomena into university classes, characterized by sets of critical exponents which encode 
information about the interface morphology and dynamics. In consequence, two different 
processes lying in the same universality class are supposed to have the same underlying 
physical mechanism driving the growth. This methodology has been successfully applied to 
a broad range of situations, but however, there are a number of restrictions that limit the 
traditional form of scaling analysis. One of them is the assumption that the interface can be 
described from a Cartesian reference frame, standard for planar surfaces. Another common 
assumption is the time independence of the substrate size. Despite the usefulness of the 
Cartesian representation in many cases, there are some growth profiles that can not be de- 
scribed according to it. Physical settings such as fluid flow in porous media grain-grain 
displacement in Hele-Shaw cells j^j, fracture dynamics j^], adatom and vacancy islands on 
crystal surfaces [4j], and atomic ledges bordering crystalline facets present interfaces 

that violate the hypothesis of the Cartesian representation. Biological systems are also char- 
acterized by an approximate spherical symmetry: bacterial colonies Q], fungi Q], epithelial 



cells [9j , and cauliflowers 



10( develop rough surfaces which are not describable from a planar 



reference frame. Also, different contexts like the technological liquid composite molding 



111, 



geological processes as stromatolite morphogenesis 12j, and chemical structures [13| provide 
examples of interfaces that either become larger as time evolves or have a curved geometry, 
revealing the broad presence of this phenomenon in the natural world. 

The theoretical study of non-equilibrium radial growth probably started with the sem- 



inal work of Eden 



14|, but the use of stochastic growth equations appeared only more 



recently, with the introduction of the Kardar-Parisi-Zhang (KPZ) equation in reparametriza- 
tion invariant form [151 ] . Subsequent works were devoted to the radial (1 + Id) KPZ equa- 



tion 



13, 



18], the radial (1 + Id) and spherical (2 + Id) Mullins-Herring (MH) equa- 
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tion 



19 



20 L and the general reparametrization invariant formulation of stochastic growth 



equations 



2l|. An analytical approach to these equations [18( showed that for short spatial 



scales and time intervals the dynamics of radial interfaces was equivalent to that of the 
planar case; however, long time intervals yielded a different output. In a recent work these 
results have been expanded showing that for fast surface growth dilution acts as a dominant 
mechanism in the large scale 22j, explaining the discrepancies in this limit. Dilution refers 
to the fact that the already deposited material becomes distributed in a larger area, and it 
thus gets diluted inside the material being deposited at the time. Herein we will limit our- 
selves to large spatial scale properties of solutions to equations whose planar counterparts 
are linear, and we will generalize our previous results [23] . The numerical study of the Eden 
model in [3] suggests that dilution is not taking place in this particular model; using this 
result as a benchmark we will phenomenologically neglect dilution in the equations under 
consideration. 

The rest of the paper is organized as follows: in section [Til we study growing radial 
interfaces and section III II is concerned with spherical interfaces; in both cases we consider 
a broader range of relaxation mechanisms than in the previous work [23]. Section HVl is 
devoted to different geometries of physical interest, which are introduced in this work for 
the sake of completeness, and finally, in section [V] we draw our main conclusions and discuss 
possible lines for future research. 



II. RADIAL GEOMETRY 



The equation of growth of a general curved surface reads [21 ] 



d t f(s, t) = h(s, t)r[f(s, £)] + $(s, t), (1) 

where the d+1 dimensional surface vector r(s, t) = {r a (s, runs over the surface as 

s = {s l }f =l varies in a parameter space. In this equation n stands for the unitary vector 
normal to the surface at r, T contains a deterministic growth mechanism that causes growth 
along the normal n to the surface, and $ is a random force acting on the surface. The 
stochastic term acts on the normal direction too $ = hrj, while the noise rj is assumed to be 
a Gaussian variable with zero mean and correlation given by 

( V (s, t) v (s\ f )> = n a (s, t)n p (s, t)eT^ 'g^ 2 5(s - s')5(t - t% (2) 



where g denotes the metric tensor determinant and the Einstein summation convention 
has been adopted. The tensor T a/3 specifies the geometric properties of the noise. In our 
case we will only consider isotropic noises, i. e. T af3 = 5 af3 , implying n a (s, i)ng(s, t)5 af3 = 1. 
Anisotropic noises appear, for instance, when considering deposition from a collimated beam 
of particles [211 ]. 

For simplicity, let us start considering radial geometry, characterized by the position 
vector 

r = (r(0,t)cos(0),r(0,t)sin(0)) . (3) 

Most of the two-dimensional examples mentioned at the introduction belong to this class: 
adatom and vacancy islands on crystal surfaces, bacterial colonies, fungi, ... . The shape 



of radial interfaces was illustrated by means of numerical simulations in 16] and 13]; their 
deterministic evolution is analogous to that of an elliptic interface of null eccentricity, see Sec. 
IIVI and Fig. [31 The effect of stochasticity in this 1 + Id system was estimated logarithmic 
231 ] . This result was obtained by means of a small noise expansion of the radial 



on average 



random deposition equation at first order in the fluctuation intensity. Of course, small noise 
expansions are appropriate only for limited times, as the unbounded nature of the Gaussian 
noise makes them invalid for arbitrary long times, for which the probability of rare events 
is not negligible. However, in our particular cases, this type of expansions yields reasonable 
results, as we will show. We will perform now a more detailed analysis of this phenomenon 
using Edwards- Wilkinson (EW) dynamics. The radial EW equation is the following Ito 



equation 



3, 



13, 



23] 



d t r = -,d A e r -- + F + — v (9, t), 

f \/V 

where r}(9,t) is a zero mean Gaussian noise whose correlation is given by 

( V (8,t) v (6>,t')) = e5(6-6')5(t-t'). 



(4) 



(5) 



Using now the second order small noise expansion 24] 



r{6,t) = R{t) + V~epi{0,t) + ep 2 {e,t), 



(6) 
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yields 



d tPi = J^Vi + Pi) + ^(M, (8) 
dtP2 = |r (fe - |p*i + P2) - 2§7^(M, (9) 



where 



(Z(6,t)ae',t'))=5(6-6')5(t-t'). (10) 

One can see that the fluctuations present in the system in the large scale come from the p\ 
dynamics. The reason is that the noise sources in the p2 equation are subdominant. This can 
be seen by deriving the equation for the correlation (p2(9,t)p2(0',t')). The most stochastic 
term appearing in the resulting equation is proportional to the four points function 

( Pl (6 1 t) Pl (6' 1 t')ae,me\t')} = ( Pl (6 1 t) Pl (6',t')} mtW,?)), (ii) 

where the equality holds because the other products of two points functions vanish due to 
the noise interpretation. This is a consequence of the Ito interpretation in Eq. (Tj0), which 
implies (rrj) = (r)(r]). Since we are decomposing the solution in the small noise expansion 
(J6]) the statistical independence among the solution and the noise holds at every order. We 



know the value of the correlator 



23| 



( Pl (9,t) Pl (9\t)) = 6 -^-P-Ht), (12) 

which increases as the logarithm of time. As the four points function is multiplied by 
R~ 3 (t) = (Ft)~ 3 , this is equivalent to have the noise term damped as in the two dimensional 



situation (see 23[ and Sec. IIIII) . making it subdominant, and the p\ random variable divided 
by \fi. The amplitude of this random variable increases as the square root of the logarithm of 
time, so this ratio vanishes in the long time limit too. The rest of the stochastic terms in the 
equation for the correlation (p2{9, t)p2{d', t')) either vanish due to the noise interpretation or 
have logarithmic amplitudes. As all these terms are multiplied by a negative power of the 
temporal variable they vanish in the long time limit, just like the terms we have considered 
here. This implies that the noise terms only contribute to create some constant prefactor, 



like in higher dimensional radial equations 23(|. This remains true for all higher orders in 



the small noise expansion, so we can conclude that the stochasticity present in the system 



comes from either the perturbation p\ or non-analytic effects. To estimate the effect of these 
exponentially infrequent events we make use of large deviation theory in appendix [Aj where 
we show evidence pointing to their irrelevance using the radial random deposition equation. 
So combining these results suggests that the stochasticity present in the system comes from 
the first order term in the small noise expansion. Of course, this result is not rigorous, 
but indicates that only a very slow process can modify it. Taking into account the explicit 
temporal dependence in the equations for the stochastic perturbations it is clear that this 
type of process will be more and more infrequent as time evolves. In the hydrodynamic limit, 
the probability of observing such a rare event will be negligible for evolution times of interest. 
This will simplify the forthcoming analyses. Note that the present situation is completely 
different to that of the KPZ or other nonlinear planar stochastic growth equations. For 
them the small noise expansion is only valid for short times, but higher order corrections in 
this expansion or a WKB analysis show the measurable effect that small fluctuations have 
on the dynamics; contrarily, these sorts of analysis have shown negligible corrections in the 
current case, as we have shown in this section and in appendix El 

Before moving to different dynamics, let us clarify some properties of the radial EW 
equation. When we talk about the irrelevance of the infrequent events we refer to the 
behavior in the hydrodynamic limit, characterized by t — > oo. For small sizes of the cluster, 
stochastic dynamics might be affected by the fixed point of Eq. (J7|), that separates the 



shrinking and the expanding tendencies 23]. In this regime, rare events might promote 



changes among these two types of evolution. In the long time limit and large spatial scale, 



however, radial EW dynamics reduces to radial random deposition 23| . and thus infrequent 



events become irrelevant, in agreement with the results reported in appendix |A] The exact 



solution of Eq. (J7j) and a study of its dynamics is reported in 23) and in Se c. IHIl and 
the analysis of the radial EW equation metastability properties is performed in [16j| and in 
appendix [0 



It is easy to generalize the EW equation to fractional orders in the planar case 25 ] 



d t h = D c A c h + ri(x,t), (13) 

where A**h denotes a Riesz derivative, and its action on the function height h(x, t) is defined 
by means of its Fourier transform 

Mh = -\kfh{k,t), (14) 



where k denotes the Fourier variable. The case ( = 2 corresponds to the EW equation. The 
object of introducing these equations is twofold. On one hand they allow understanding 
the interplay of the value of exponent (, which specifies the velocity at which correlations 
propagate, and the other characteristics of the system. On the other hand, these fractional 
operators are common models describing anomalous transport properties of complex systems 



261 ]. In this sense, one can think of fractional equations as physical generalizations of the 



classical Laplacian diffusion. It is not so straightforward however to obtain these equations 
in a general manifold. Let us recall the Beltrami-Laplace operator in one dimension 

w = > O f ) ■ (is) 

where s parameterizes the corresponding manifold. Its structure suggests us a possible 
definition of "square-root" of this operator as 

Af=^-A s f, (16) 
V 9 

where A s is some derivative with the desired properties. Due to the intrinsic nonlinearity 
of this operator, we cannot define a fractional derivative along the lines of the Riesz one, 
because it takes advantage of the linearity of the ordinary derivative in order to use the 
Fourier transform. Nevertheless, we can generalize the radial EW equation to fractional 
orders assuming r(9,t) = R(t) + y/ep(9,t), R(t) = Ft, and the perturbation obeying the 
equation 

d tP =j^Aip + -±=ae,t), (i7) 



in the interval [0, 27r], where the solution is subject to periodic boundary conditions. This 
equation, in the long time limit, coincides with the corresponding EW one when £ = 2, 
except for one term coming from the nonlinear dynamics of the curved surface, see Eq. (JS]). 
So this toy model retains part of the properties of the equation of motion for the inter- 
face, but part of the curvature effects is lost, and these are the cause of some stability 
deterioration as have been seen 13, |23|]. However, we will see that the correlation ob- 
tained with this model is the same that the one obtained from the original equation for the 
perturbation. The correlation can be computed in terms of the Fourier transformed field 



7 



p n {t) = (27T)- 1 P (9,t)exp(-in9)d9: 



(p n (t)p m (t)) = exp 



D c (\m\ 



\n\ 



exp 



(Pn(t )p m (t )) + 



8 n 



2ttF(C 



Ei 



£>c(M c + 



n\ 



Ei 





"D c (|m| f + 
. FCtf-i(l-C) 



(18) 



if C 7^ 1 (note that for ( = 2 the result is the same as for the full radial EW equation 23|); 
here n,m G Z, Ei(x) = f^ oQ t~ 1 e t dt denotes the exponential integral 271 ]. and this last 
integral has to be interpreted in terms of the Cauchy principal value when necessary. Using 
the asymptotic expansion Ei(sc) ~ ln(|a;|) when x ~ we find, when £ > 1 and in the limit 
t — > oo, that the correlation adopts the form 



(Pn{t)p m {t)) 



6„ 



2ixF 



Mt), C(9,9',t) = (p(9,t)p(9',t)) 



2ttF 



5(9-9'), 



(19) 



and so reduces to radial random deposition, this is, becomes totally uncorrelated. If ( < 1, 
and in the long time limit, we find that the correlation vanishes 



(p n {t)p m {t)) -> 0, 



(20) 



by means of the asymptotic expansion Ei(cc) ~ e x /x when x — > ±oo, and the resulting 
interface is flat. If ( = 1, what we call the ballistic model, the correlation is 



1 

2^ 



C(9,9',t) = 
sinh[2Di(r - s) 



2tt 



e 2|n[Ci( S -r)+in(0-0') rfs 



-ds 



:ln 



t 



(21) 



cosh[2Di(r-s)] - cos(9 - 9')^ 2irF~* \\6 - 6'\ F / D \ 

where r = \n(t)/F. The final expression has been obtained in the limit of long times and 
small angular differences. This result allows us to define the local dynamical exponent 
zioc = F/Di, which depends continuously on the ratio of the parameters F and D\, and it is 
thus nonuniversal. The planar ballistic model is analyzed in appendix [UJ where it is shown 
that the correlations in this case propagate indeed ballistically, this is z = 1 is the universal 
dynamical exponent. 

The radial MH equation reads 



19 



20] 



d t r = ~rM r + d e r ) + F + -4^ °- 1 1 • 



(22) 



8 



It appeared due to its possible connection to tumor growth [19| . By means of a small noise 
expansion we find at the deterministic level, compare to Eq. (j7]), R(t) = Ft, while for the 
stochastic perturbation, compare to Eq. we get 



dtp = -j^-( d eP + d oP) + -fiffr *)■ 



(23) 



The Fourier transformed field obeys the equation 

dp n K 



dt 



FH 4 



(n -n )p n + —=n n (t) 



Ft 



(24) 



that yields the mean value 

(Pn(t)} = (p n (t ))exp 



K(n 2 - n 4 ) ( 1 1 



3F 4 



t\ *3 



(25) 



showing that all the modes are stable except n = 0, ±1, which are marginal, a situation rem- 
iniscent to that of the EW equation. The equation for the correlation C n , m {t) = {p n (t)p m (t)) 



da 



n,m 



dt 



K(n 2 + m 2 — n 4 — m 4 ) 5 n , 



FH 4 



2nFt 



(26) 



can be solved to get 



C nt m(t) = exp 



K{n 2 + m 2 — ri 4 — m 4 ) ( 1 1 



Cn,m(^o) + 



Ei 



6ttF 

Kin 2 + m 2 — n 4 — m 4 
3^ 



exp 
- Ei 



3F 4 V^o t3 

K(—n 2 — m 2 + n 4 + m 4 
3FHI 

K{n 2 + m 2 — n 4 — m 4 ) 
3FH 3 



x 



(27) 



Assuming an uncorrelated initial condition and using the asymptotic equivalence Ei(x) 
In (\x\) when x ~ yields 



(Pn(t) Pm (t)) ~ feln^), <7(0, t) ~ - 0') 



2ttF 



2ttF 



(28) 



when £ — > oo. We have found that in this case again, as happened with all the models such 
that z > 1, the interface becomes uncorrelated in the long time limit. 

In order to unify the different equations let us build a generic one assuming that the drift 
comes from a potential 

r ^ m 1 SV\r(s,t)] . . 

r[f( s , t)] = - _ jLyJ , (29) 



y/g{s) Sf(s,t) 



9 



and this potential can be expanded in a power series of the surface mean curvature H 



» N N 

J i=0 i=0 



(30) 



The general contribution to the drift then reads 



r< = — 



Ki I H i+l - tAsLH 1 - 1 



3=1 



(31) 



where Xj are the eigenvalues of the matrix of the coefficients of the second fundamental form 
and express the principal curvatures of the surface. Applying this theory to the present case 
of radial growth we get an equation generalizing both EW and MH 

' dlr r 



d t r = K Mr - - 



1 f) 4 r \ 1 



(32) 



where we have linearized the resulting equation in the different derivatives of r about zero 
(limit of small variations of the radius with respect to the angle) and have truncated the 
expansion at order two; note that this procedure left us with the EW (zeroth order, propor- 
tional to Kq) and MH-like (second order, proportional to K 2 ) terms, while the first order 
term vanishes identically (a different situation arises in higher dimensions, see section fTTTT) . 
This equation is particularly interesting for its stability properties. The usual small noise 
expansion leaves us with a system of equations for the different Fourier modes, whose mean 
value is 

(t - t )[3F 2 K (m 2 - l)tHl + K 2 (3 + 7m 4 )(t 2 + tt + t 2 )} 



(Pm(t)) = exp 



3FHH 3 



(Pm(to)} 



(33) 

revealing that the modes m = ±1,±2, ... are stable and thus improving the stability of 
both EW and MH equations, for which the modes m = ±1 were only marginal 23(|. The 
stochastic properties of this equation are also attractive from a metastability point of view, 
see appendix (Bj what suggests it is a reliable model for describing radial growth processes. 



III. SPHERICAL GEOMETRY 



In this section we concentrate on 
appear when considering tumor growth 



spherical (2 + Id) growth models. These models can 
or other biological structures, or diffusing voids 



3 
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FIG. 1: Deterministic growth of the spherical cluster at zeroth order in the small noise expansion, 
see Eq. (|36p in the main text, with parameter values F = 1, Kq = 1/2, and to = 1 + 10 -4 (red 
lower line). The linear law t + 1 + 10 -4 is plotted for comparison (green upper line). 



inside a solid [28]. We will analyze the basic properties of the simplest equations that arise 
in this geometry using the position vector 



f = (r(0, 0, t) cos(0) sin(0), r(0, <f>, t) sin(0) sin(0), r(0, 0, t) cos(0)). 



The spherical EW equation reads 



dlr 



F + 



ri(0,<i>,t). 



(34) 



(35) 



r 2 tan(0) ' r 2 ' r 2 sin 2 (0) 

Note that its drift is proportional to the mean curvature of the interface T = K H, once 
the small gradient expansion has been performed. As usual, we will decompose the solution 
into a deterministic radial part and a small stochastic perturbation r(t) = R(t) + ep(8,t). 
The deterministic part obeys an equation that is analogous to the one obeyed by its 1 + Id 



counterpart R = F — 2K R x , see Eq. ([7]), and that can be solved to yield [23] 



R(t) 



2K n 



1 + Wo 



FH 



1 exp 



FH 



(36) 



2K J r \2K 

where the initial condition was assumed to be R(to) = Fto, and it behaves asymptotically 
in time as R(t) ~ Ft when F 2 t > 2K , see Fig. [TJ if this inequality is reversed then the 
solution shrinks till it collapses in finite time 23j , see Fig. [21 Here Wo denotes the principal 
branch of the Lambert omega function 29(. The random perturbation obeys the equation 



dtp 



K 
FH 2 



9dP +%p+ 91 



tan(0) 



sin 2 (0) 
11 



+ 2p 



+ 



Fty^ue) 



77(0, M- 



(37) 



Now we assume that the solution can be decomposed as a linear combination of the spherical 
harmonics 

oo I 
1=0 m=-l 

The linearity of Eq. fl37|) allows us to write the evolution equation for the different modes 

«-£ P -,<« + 1,^*0. (30) 

where the spherical noise has been expanded in spherical harmonics as well 

rj(9,<f>,t) = 

1=0 m=-l 



oo I 



E E ^rw r 



(40) 



Vsin^ 

and 77™ (t) is a zero mean Gaussian distributed random process, whose correlation reads 

Note that the presence of the Condon-Shortley phase in the correlation makes the noise 
change from real to imaginary and back as we vary m. We can solve Eq. fl39l) for its mean 
value 



(pT(t)) = exp 



Ko(l 2 + l-2)(t-to) 
FHta 



(P?(t )) 



(42) 



which is independent of the value of m except for the initial condition, and shows the linear 
stability of the spherical symmetric phase for the values I > 1. Perturbations with I = are 
unstable, while I = 1 characterizes the marginal case, a situation reminiscent to that of 1 + 1 
dimensional radial growth [23]. Its correlation function C^; m (t) = (p/ ra (^)P™'(^)) obeys the 
differential equation 

. — l) m ^m,-m'^,/' 



dc^r _ k [a - i(i + 1) - + , 



dt 



FH 2 



F 2 t 2 



(43) 



that can be solved to yield 

i,v (*) = C l,i> (*o) exp 
1 — exp 



K [l{l + l) + l'(l' + l)-4){t-t ) 



+ 



FHt c 

K [l(l+l)+l'(l'+l)-4](t-t ) ~ 
F 2 tt 



(44) 



K [l(l + l) + l'(l' + l)-4] 
The correlation is bounded as t — > oo, and it is a function of to- If we now take the limit 
to — ► oo we find 

lim \imC-r'(t) = 0. (45) 

to— >oo t— >oo ' 



12 




FIG. 2: Deterministic collapse of the spherical cluster at zeroth order in the small noise expansion, 
see Eq. (|36p in the main text, with parameter values F = 1, K$ = 1/2, and io = 1 — 10 -4 . 

This means that fluctuations average out in the limit we are considering. The physical 
reason is that the surface area continuously grows while the parameter space stays constant. 
However, the surface might accumulate some average roughness produced by perturbations 
in earlier stages of growth. This finite size effect reflects the memory of the growth process 
with respect to the time the perturbation was set in. This is a consequence of the existence 
of an absolute origin of time in this system: the instant characterized by a cluster with 
zero radius. As the temporal and spatial scales of this system are intimately related, the 
hydrodynamic limit is characterized by an infinite lapse of time from the absolute temporal 
origin. Interestingly, perturbations happening at the first stages still remain in the interface 
in this limit; however, perturbations starting later have no effect on the average dynamics. 
We have shown elsewhere that this effect is a consequence of the absence of dilution and a 



fast radius growth 



22]. 



If we perform an expansion from a potential in a power series of the mean curvature, just 
like at the end of section [Til we wm nn d again the spherical EW equation at zeroth order, 
but in this case the first order does not vanish. If we use this first order as the drift of a full 
stochastic growth equation we obtain what we call the intrinsically spherical (IS) equation 



d t r = K x 



~&kr , dlr d e r 



+ 



F+ -L= q(M,*), (46) 
ry sm(0) 



r 3 r 3 sin 2 (6 l ) r 3 tan(#) r 2 

that has no 1 + 1 dimensional analog. We proceed to analyze this equation as usual. The 
solution of the spherical symmetric equation R = F — K\R~ 2 can be obtained in implicit 



13 



form 



R(t) 




arctanh 




R(t) 



Ft- 




arctanh 



F_^% 



(47) 



The solution grows unboundedly for initial conditions such that F 3//2 t > K]J 2 , and it adopts 
the linear growth R ~ Ft in the long time limit; when this inequality is reversed then the 
solution shrinks till it collapses in finite time. The equation for the stochastic perturbation 
reads this time 



d t p 



FH 3 



d e p , a2_ , ®\p 



tan(0) 



+ d 2 eP + 



sin 2 (0) 



+ 2p 



+ 



1 



Fty^m(9) 



ri(0,<i>,t). 



The spherical harmonics decomposition yields for the mean value of the perturbation 



Of(t)> = exp 



K 1 (l 2 + l-2)(t 2 -t 2 c 
2FHH 2 



(48) 



(49) 



what reveals that the mode / = is unstable, the modes with I = 1 are marginal, and 
the rest of the modes is stable, exactly the same as for the spherical Edwards- Wilkinson 
equation. For the correlation function we get 

\ffi[Z(Z + l) + /'(/' + !)- 4] 



C U> W = eX P 



2F 3 

+( — l) m ^m,-m'^y 



it- 2 -to 2 ) 



FK^A- 1(1 + 1) + 



X 



exp 



K x [l(l + \) + V(V + l)-4\ 2 a 



2F 3 



^[4 - Z(Z + 1) -/'(/' + !)] ^ 
2F 3 f 



^[4 — Z(Z + 1) — Z^CZ^ + 
2F 3 



-V 



where T>(x) = e x2 j Q x e y2 dy is the Dawson integral [27|. Using the fact that T>(x > ~ 



(50) 



x when 



x ~ 0, we can deduce that the correlation is bounded as t — > oo, and we also recover the 
same result, Eq.(j4"5|). as in the last case. 

The spherical Mullins- Herring equation reads [3, Q 

d t r = --| {[2 + sin~ 2 (^)] tan -1 ^)^ - tan- 2 (#)<9 e 2 r + 2 tan" 1 ^)^ + d 4 e r 
-2 sin~ 2 (#) tan -1 ^)^^ + 2 sin- 2 (fl)<9j<9 2 r + sm- 4 (#)[4<9jr + <9jr]} 

+F+ J_ q(g,6*)- (51) 
rysm(&J 
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It appears when surface diffusion occurs to minimize the surface area, yielding the drift 
term T = — K^blH ■ In this case, our usual division of the solution gives the very simple 
expression for the deterministic radial part R(t) = Ft, while the stochastic perturbation can 
be expressed as an infinite series of spherical harmonics. For the different modes we get the 
equation 



dt FH 4 
and so for the mean value we obtain 



K " [z(/ + i)(z 2 + / + 2)K l + rr(t) 



Ft ' 



(p?(t)) = exp 



K 2 [l(l + l)(P + l + 2))._ 



3F 4 



(pT(to)) 



(53) 



what reveals that all the modes are stable, but I = that is marginal. For the correlation 

we get 



C^r'(t) = exp (J^[l(l + 1)(Z 2 + / + 2) + /'(/' + 1)(/' 2 + V + 2)(r 3 - V 3 )]) 



x 



err + 



3F 2 



exp 



3FH 3 



[1(1 + l)(l 2 + l + 2) + l'(V + 1)(/' 2 + l' + 2)))x 



6 m ,-m>5i,i' |*~V-2/3 {j^W + 1)(/ 2 + Z + 2) + l\V + 1)(/' 2 + + 2)]^ 

"*o V-2/8 (^PC + 1)(/ 2 + / + 2) + /'(/' + 1)(/' 2 + 1' + 2)]j } : 



(54) 



where (p n (x) = y n e~ xy dy is the Misra function 36| . Using the expansion yj_2/3(x) ~ 
r(l/3)x -1 / 3 (here T(x) denotes the gamma function and T(l/3) ~ 2.68) when x ~ we 
see that the correlation is bounded in the infinite time limit and that we recover the same 
result, Eq. (j4"5]) . as in the other two spherical cases. 



IV. OTHER GEOMETRIES 



It is also interesting, both at the theoretical and applied levels, to consider geometries 
that do not preserve the spherical symmetry. The simplest case might be elliptic geometry, 
see Fig. [31 An example of elliptic interface is provided by epithelial cell populations on 
stretched elastic substrates We will consider the position vector 



r = (r(8,t) cos(6),ar(6,t) sin(0)), 



(55) 



where a > 1 is some parameter specifying the ellipse eccentricity. Following the techniques 



developed in 



2l|, we can derive the EW equation in elliptic geometry. However, as could 
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FIG. 3: Sketch of a deterministically growing elliptic interface of constant eccentricity. Red (inner) 
line, short times; green (middle) line, intermediate times; blue (outer) line, long times. Stochasticity 
will manifest itself as small Gaussian fluctuations about the deterministic shape in the small noise 
approximation . 

be expected, the resulting equation is not spatially homogeneous, as it explicitly depends 
on the angle 9. As an example, we show the expression of the metric tensor (note that, due 
to the low dimensionality, the metric tensor is a scalar in this case), keeping only the terms 
that are linear in the derivatives of the radius 

g = a 2 r 2 cos 2 (9) + r 2 sin 2 (9) + (a 2 - l)rr e sin(20). (56) 

For the sake of simplicity, we will derive two EW equations, one that locally describes the 
points that are equidistant to the ellipse foci, which is 

er \ r z r J ^ar 
and the other one for the points that lie the closest to one of the foci 

/ Ffir 1 \ F 1 

*-^--) + - + ^(M. (58) 

The advantage of these concrete points is that they are the only ones affected by isotropic 
growth. As can be seen, the growth is qualitatively identical to that of the radial cluster. 
At the quantitative level, one can see that the neighborhood of the minor axis is affected by 
a reduced diffusion, while the neighborhood of the major axis is affected by reduced growth 
and fluctuations. Let us now consider the random deposition model in ellipsoidal geometry 

f = (ar(9, 0, t) cos(0) sin(0), r(9, 0, t) sm(9) sin(0), r(9, 0, t) cos(0)); (59) 
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FIG. 4: Sketch of a deterministically growing hyperbolic interface of y/2 eccentricity. Red (inner) 
line, short times; green (middle) line, intermediate times; blue (outer) line, long times. Stochasticity 
will manifest itself as small Gaussian fluctuations about the deterministic shape in the small noise 
approximation . 

we will concentrate on the dynamics on the 8 = tt/2 plane. The ellipse lying on this plane is 
similar to the one just described above but the random deposition model reads in this case 

d t r = - + -aO,<f>,t), (60) 
a ar 

near the major axis, and 

d t r = F + -^£(6,<t>,t), (61) 
ry a 

near the minor axis, where £(#, 0, t) is a zero mean Gaussian noise whose correlation is given 
by 

(£(0, 0, t)£(6K, 0', t')} = e5(6 - 9')8{4> - <f>')6(t - t'). (62) 

Like in the 1 + Id case the neighborhood of the major axis is affected by reduced growth 
rate and fluctuations intensity, while the neighborhood of the minor axis undergoes a smaller 
reduction of the fluctuations intensity but its growth rate remains unchanged. As in the 
spherical case, the strong negative power of the radius in the noise term will cause a constant 
average roughness that depends on the perturbation initial time. 
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The opposite case to elliptic geometry is hyperbolic geometry, see Fig. @J Thisgeometry 
could perhaps be useful in describing atomic ledges bordering crystalline facets [5j. In this 
case, different layers in the solid corner would correspond to the different hyperbolas. If we 
define the position vector as 

r= (p(x,t)sinh(x),p(x,*)cosh(x)), (63) 

then p > is the "radius" of the hyperbola upper branch, to which we will limit our 
study, and \ £ K. is the hyperbolic angle. This geometry is again characterized by a 
nonhomogeneous growth equation, so we will focus on the hyperbola vertex in order to 
obtain an isotropic result. In this case the resulting EW equation reads 



%P = D[^ + l)+F + -^t(x>t)- ( 64 ) 



Performing as usual the expansion p{x,t) — Tlif) + ecr (x ? ^) ; where e is the small noise 
amplitude, we get the zeroth order equation 



dTZ _ D 

lu = n' < 65) 



that can be solved to yield 



K(i) ^J 1 + nu(-£±^exp 



D + F(JZ + Ft) 
D 



(66) 



where W-i is the corresponding branch of the Lambert omega function [221], and IZq = 71(0). 
One can read from this equation that the growth is faster than linear, but it evolves towards 
the linear law 7Z(t) ~ Ft in the long time limit, see Fig. [5j Furthermore, there are no fixed 
points in the dynamics, and the system can only expand, and not shrink like in the radial 
case. The equation for the perturbation is 

a <° =D {w)-m) + WY (67) 

where as usual r\ is a zero mean Gaussian noise which correlation is given by 

(v(x,t)v(x',t')) = e5(x-x')S(t-t f ). (68) 

Now one sees that, contrary to the radial case, all the terms in the drift of the equation 
for the perturbation are stabilizing, and even a homogeneous perturbation decreases in time 
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FIG. 5: Deterministic growth of the hyperbolic interface at zeroth order in the small noise expan- 
sion, see Eq. (|66p in the main text, with parameter values F = 1, D = 1, and TZq = 10~ 4 (red 
upper line). The linear law t + 10 -4 is plotted for comparison (green lower line). 

due to their action. To study the higher dimensional growth properties of a hyperbolic 
surface let us now consider a double sheeted revolution hyperboloid. The upper sheet may 
be parameterized with the vector 

r = (p(x, 0, t) sinh(x) cos(0), p(x, 0, t) sinh(x) sin(0), p(x, 0, t) cosh(x)), (69) 
and the random deposition model, next to the vertex, reads in this case 

d t p = F+—^=axA,t), (70) 
PV\X\ 



what reveals that the hyperbolic surface develops a constant average roughness in this case 
as well, just like the spherical and ellipsoidal geometries. Note that the vertex of the hyper- 
boloid, parameterized as x = 0, is a singular point characterized by an infinite amplitude 
of the fluctuations. This paradox is solved by taking into account nonlinear terms in the 
derivatives of the field 

dtp = F + 1 gfr, 0, t). (71) 

(/xV + p 2 pI 

Finally, another interesting geometry is that of a growing parabolic surface, see Fig. El 
Stromatolite morphogenesis might be related to parabolic geometry, as has been pinpointed 



in recent works 



30 



3jj. Using the parameterization 

f = (2x$(x,t),(l -x 2 )$(x,t)), (72) 
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where x G M and $ > 0, we obtain the EW equation 

a <*= D @-2¥) +F+ vb« M >< (73) 

for the dynamics close to the vertex of the parabola, x = 0, the only point where isotropy 
holds. The noise £ is Gaussian and zero centered, and its correlation is given by 

(C(x, t)f 0) = - x')5(t - t'). (74) 

One can see that this equation describes a dynamics which is identical to the spherical one, 
up to some numerical adjustment. The paraboloidal geometry can be characterized by the 
position vector 

f = (2x$(x, 0, t) cos(0), 2x$(x, 0, t) sin(0), (1 - x 2 )$(x, 0, £)), (75) 

that rends a random deposition model 

d t <S> = F+—± 7 =Z{x,<j>,t), (76) 

where the zero centered Gaussian noise is again delta correlated and has some intensity e. 
As in the hyperbolic case, the noise produces an average constant roughness depending on 
the perturbation initial time, and the vertex is again a singular point. Taking into account 
the lowest order nonlinearity in the vertex the resulting equation is 

<9 t $ = F + 1 =e(x, 0, t). (77) 



As can be seen, parabolic interfacial dynamics lies between the elliptic and hyperbolic cases. 
This is rather natural, as a parabola is the conical curve that represents the marginal situ- 
ation between the ellipse and the hyperbola. 



V. CONCLUSIONS 

The dynamics of radially symmetric surfaces has been studied in sections [TT] and IIHI In 
Sec. IIHI we have calculated for the first time the explicit correlations of the spherical models, 
and we have shown that the long time average roughness becomes constant in polar coor- 
dinates, but at the same time it develops memory with respect to the initial conditions. In 
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FIG. 6: Sketch of a deterministically growing parabolic interface. Red (lower) line, short times; 
green (middle) line, intermediate times; blue (upper) line, long times. Stochasticity will manifest 
itself as small Gaussian fluctuations about the deterministic shape in the small noise approximation. 

Sec. [IT] it is shown also for the first time that sub-ballistic diffusion radial models reduce to 
radial random deposition in the long time limit and for large angular scales. Radial random 
deposition is characterized by an uncorrelated interface of a logarithmic temporal ampli- 
tude, as was shown in [23]. The model with ballistic correlations propagation, introduced 
here in radial geometry, shows a richer phenomenology. Despite the logarithmic temporal 
dependence, the correlation function is markedly different. In fact, the analysis of this func- 
tion for long times and short angular scales reveals the appearance of a local dynamical 
exponent zi oc , which is equal to the ratio of the constant rate of radius growth F and the 
ballistic diffusion constant D\, and it is thus nonuniversal. This reflects the fact that the 
microscopic details of the growth process influence the surface scaling at the local level. We 
have also shown for the first time that our results in radially symmetric geometries seem 
to be metric independent. We have been able to reproduce them in elliptic, hyperbolic, 
and parabolic geometries, which undergo particular nonlinear effects, but share the same 
stochastic dynamics. The differences among the classical results on surface growth and the 
current ones are related to the substrate growth properties. The classical setting assumes 
a constant d dimensional substrate size, and growth is restricted to the (d + l)th spatial 
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dimension. In the parameterizations we have considered here, growth happens equally in 
every spatial dimension, a fact that has been ensured by keeping constant the eccentricity 
of the different underlying conical structures. Indeed, in Sec. IIVI we have kept the eccen- 
tricity constant explicitly by means of the introduction of the parameter a. In this same 
section, hyperbolic geometry has been implicitly supposed of \/2 eccentricity, but assum- 
ing an arbitrary constant eccentricity would yield a straightforward modification of these 
results, just along the lines of the elliptic case, and it would rend the stochastic dynamics 
unchanged. A different situation might arise if the eccentricity changes along the evolution. 
For instance, consider an ellipsoid growing in all spatial dimensions at different rates. Its 
shape will change in time, and this might affect the surface dynamics. In this case, different 
spatio-temporal scales will enter in competition, and the resulting growth dynamics might 
combine features of the different scenarios studied so far and perhaps new ones. In any case, 
due to the complex nature of these growth regimes, it is necessary to build some analytical 
progress before drawing any conclusions about them. 

A necessary step in the following is the study of the relationship among stochastic growth 
equations and discrete models. On the applied side, important technological processes such 



as liquid composite molding ll( are driven by growing radial interfaces. The theoretical 
study of this sort of systems is often based on detailed numerical simulations, which could 
be simplified using eikonal descriptions like the stochastic growth equations presented here. 
The duality expansion versus shrinking present in the radial and spherical EW equations 
and IS equation might be related to some of the observed phenomena. Interface shrinking is 
present in adatom islands, which can disappear if they are composed of a subcritical number 



of particles 32] . The comparison with numerical simulations and experiments will facilitate 



a deeper understanding of the dynamics of curved surfaces. 
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APPENDIX A: ESTIMATING THE EFFECT OF EXPONENTIALLY INFRE- 
QUENT EVENTS 



In this appendix we will estimate the effect of rare events on the evolution of the radial 
random deposition process using large deviation theory. This process is described by the 
radial EW equation after setting D = 



23j . The equation of motion in the one dimensional 

case reads 

d t r = F + r- l/2 r](e,t), (Al) 
that is actually a stochastic ordinary differential equation. The noise correlation is given by 

( V (6,t) v (d',t')) = e5(6-d>)5(t-t>). (A2) 

We will assume the following form of the probability distribution (which is in fact a WKB 
ansatz) 

P(r,t) = exp[Q(r,t)/e], (A3) 
and in the limit of vanishing noise intensity e->0we find 

d t $ = -Fd r $ + (A4) 

The equation for its derivative \1/ = <9 r $ can be solved along characteristics 

Tt = (A5 » 
Tt =F ~7- < A6 > 

The analysis of the characteristic equations show the temporal behavior of the radius 

r (t) «F* + fen(Ft), (A7) 
b 

in the long time limit, where ipo = ^(0) measures the initial size of the rare fluctuation, and 
we have assumed that Fr(0) 3> l^ol- This result shows that the effect of these exponentially 
infrequent events is negligible in the long time limit. Redoing this same analysis for d > 1 
one finds that the logarithmic correction is substituted by a constant proportional to ipo, 
what shows that the effect of rare fluctuations is even weaker in this case. This analysis 
totally agrees with the results obtained from the small noise expansions in Sec. [IT] and IIHI 
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V 




FIG. 7: Deterministic (green lower line) and effective stochastic potential (red upper line) for the 
mean- field model of the radial EW equation. The parameter values are F = 1, D = 1, and e = 0.4. 

APPENDIX B: METASTABILITY IN THE RADIAL EDWARDS- WILKINSON 
EQUATION 



In Ref. 1_6| the authors investigate the behavior of the interface driven by the radial EW 



equation. To this end, they use the Ito stochastic differential equation 



%=F -- + -!=&)> ( B1 ) 
at r Jr 



(tmt'))=e6(t-t% (B2) 

as a mean field model, and they showed that it represents well the dynamics of a lattice 
simulation interface. The deterministic fixed point is given by r c = D/F, that separates 
the shrinking and growing phases. However, the presence of the noise term can alter this 
scenario. We can interpret Eq. ( IBlj) as the overdamped description of a particle moving on 
the singular potential V(r) = Dln(r) — Fr, that possesses two wells of infinite depth, the 
origin and the positive infinity, and it is subject to state-dependent noise. We can calculate 
the mean time that it takes this particle to leave the well at the origin, and to roll down to 
infinity after overcoming the maximum at r c . It is given by the solution to the well known 



equation 



24) 



—T" + ( F — — J T = — 1, (B3) 
2r \ r J 

where T = T(r) is the mean time at which a particle starting at r, < r < r c , reaches 
the position r c , and thus, the mean first passage time can be defined as 2T. The correct 
boundary condition at r c is an absorbing one T(r c ) = 0. At the origin we will assume we 
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have a reflecting boundary T r (0) = 0, a fact that we will explain in more detail below. 
The solution of this equation subject to the mentioned boundary conditions can be readily 
computed, but for simplicity we will only state the result for an initial condition at the 
origin T = lim r ^ Q T(r) 



F 2 1 
—T = 1 + -erf (a;) 



7rerfi(x) — ^-!L e x2 

x 



- 2 T 2 (l,l;3/2,2;-xV, (B4) 



where x = D/y/eF, erf (2) = 2tt~ 1 / 2 J* e~ y dy is the error function 27], erfi(z) = —ieri(iz) 
is the imaginary error function, and 2-^2 a 2] W, ^2; z ) is the corresponding generalized 
hypergeometric function {33)]. Of course, if we solve the passage problem for an initial 
condition such that r c < r < +00, the mean first passage time is always divergent. 

Despite its exactness, this formula is very complex, and a further clarification of the 
passage process is in order. To accomplish this, we can use the Stratonovich version of 
Eq. (JET]) 

dr D e 1 , . . . 

T t =F -7 + ^ + v?° m (B5) 

and by changing variables u = r 3//2 we obtain the additive noise equation 

~ = Fu^ - Du-W + V 1 + m = -\V u {u) + m, (B6) 

where the new potential is defined as 

^V u (u) = -6Du 2/3 + 3Tm 4/3 + eln(ii). (B7) 

The evolution of the probability functional associated to the random variable u(t) can be 
cast in the form of a path integral [34| which action reads 

2 f°° 

SHt)} = - {u(t) + V u [u(t)]} 2 dt (B8) 
J J —00 

Its first variation SS/Su = gives the equation for the classical trajectory 

u = V u Vl (B9) 

that admits a first integral of motion, and using the fact that the particle stops at either 
the minimum or the maximum of the potential, we conclude that the integration constant 
is zero. We are thus led to the equation 

u 2 = V>) 2 . (BIO) 
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Taking the square root we find 

l^L = F^/s _ Dw -i/3 + £ -i (Bn) 
3 dt 4 v ; 

where we have selected the negative square root. The positive square root is the instanton 



3J. Re- 



solution that indicates the optimal path for escaping from the potential well 
changing variables to recover the description in terms of r we find 

— = F- — + —, (B12) 
dt r 4r 2 ' 



which is nothing but the Stratonovich equation (1B5I) with the noise term suppressed 35]. 
The mean field equation (1B1 2j) represents motion on the effective potential 



V e = Dln(r)-Fr + -^, (B13) 
4r 

which is positively divergent at the origin, justifying the choice of a reflecting boundary at the 
origin earlier in Eq. (1B3[) . This effective potential is a regularized version of the corresponding 
one in Eq. ( 1B1I) . and it has a finite depth minimum at r m j n = (2F)~ 1 (D — \J D 2 — eF), and 



a finite height maximum at r max = (2F)^ 1 (D + v 1 D 2 — eF), provided D 2 > eF. Otherwise 
it is monotonically decreasing in r, implying the loss of met ast ability. So we see that the 
simple picture of bistability drawn by the small noise expansion can be modified due to 
the divergent amplitude of the fluctuations at the origin. However, the behavior for long r 
remains the same. A comparison between both potentials can be found in Fig. [71 

For the bistable spherical models EW and IS we can derive the same kind of mean field 
equations 

(Bi4) 

r 

h k(t), (B15) 

respectively, and the corresponding semiclassical description 

dr 2Ko e ,„ 

*= F "T 1 + ^ (B16) 

dr „ K\ e 

Tt= F -^ + ^- < B17 > 

In both cases we find the same regularization as for the radial EW equation, because we 
know, by applying Descartes rule, that in the stationary regime these last two equations 
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dr 




2K 




= F — 


~db 




r 


dr 




Ki 




= F — 









have either none or two roots. The exact values of the maximum and the minimum of the 
potential, and the threshold condition for the existence of bistability are again within reach, 
since any cubic polynomial equation can be solved in terms of radicals. However, as usually 
happens, the solutions are cumbersome and not particularly illuminating, and so we will 
omit the exact expressions here. Of particular interest is Eq. (I32p . which mean field model 
reads 




(B18) 



In this case, suppressing the noise term either in the Ito or the Stratonovich interpretation 
yields the same qualitative behavior of the effective potential. This shows that the determin- 
istic evolution of this equation is more robust to stochastic perturbations, and suggests that 
the type of equations that can be obtained from series like (1301) might be a good starting 
point for modelling the growth of complex interfaces with general geometric symmetries. 



APPENDIX C: THE PLANAR BALLISTIC MODEL 

The ballistic model in planar geometry reads 

d t h = D x Ah + £{x,t), (CI) 

and its Fourier transform is 

d t h = -Di\k\h + h(k,t), (C2) 

that can be solved to yield 

/oo 
e Dl ^ s ^i(k, s)ds. (C3) 
-oo 

The correlation in Fourier space reads 

(h(k, t)h(q, t)) = e J e Dl ^ + ^ s ^5(k + q)ds, (C4) 

and in real space 

poo /*oo pt 

(h{x,t)h(y,t)) = e / dk dq ds[e^ kx+qy) e Dim+ms - t] 5(k + q)] = 



oo J — oo J 

2+2 



2Di \ \x — y\ 2 J Di \\x — y\J ' 
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when t — > oo. As a final note, let us remark that this simple analysis has been made in 
order to obtain a comparison with the results in section [Til A more detailed analysis can be 
found in [251 ]. 
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